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ABSTRACT 

We have observed a bright flare of Sgr A* in the near infrared with the adaptive optics assisted 
integral field spectrometer SINFONI^. Within the uncertainties, the observed spectrum is featureless 
and can be described by a power law. Our data suggest that the spectral index is correlated with the 
instantaneous flux and that both quantities experience signiflcant changes within less than one hour. 
We argue that the near infrared flares from Sgr A* are due to synchrotron emission of transiently 
heated electrons, the emission being affected by orbital dynamics and synchrotron cooling, both acting 
on timescales of w 20 minutes. 

Subject headings: blackhole physics — Galaxy: center — infrared: stars — techniques: spectroscopic 



1. INTRODUCTION 



The detection of stel lar orbits (iSchodel et alJ 120021: 
lEisenhauer et l2(M IGhez et aljl 



close 

to Sgr A* has proven that the Galactic Center (GC) 
hosts a massive black hole (MBH) with a mass of 
(3.6 ± 0.2) X 10*^ Mq. Sgr A* appears rather dim in 
all wav elengths, which is explained by accretion flo\y 
models fNaravan et al."' l995HQuataert &: Naravanlll999t 
iMeha et al...2000. ,2001.) . In the near infrared (NIR) it 
was detected after diffraction limited obser vations at the 
8-m class telescopes had become possible (i Genzel et al.l 
[2003; Chez et al. 2004). Usually the emission is not de- 
tectable. However, every few hours Sgr A* flares in the 
NIR, reaching up to K k, 15 mag. A first fl are spec- 
trum was obtained bv lEisenhauer et alJ l|2005D . showing 
a featureless, red spectrum [uSi, ~ with /? « —3). 

2. OBSERVATIONS AND DATA REDUCTION 

We observed the GC on 2 005 June 18 from 2 :40 to 7:15 
UT w ith SINFONI ( Eisen hauer et al. 2003; Bo nnet et ail 
l2004|) . an adaptive optics (AO) assisted integral field 
spectrometer which is mounted at the Cassegrain focus of 
ESO-VLT Yepun (UT4). The field of view was 0.8" xO.8" 
for individual exposures, mapped onto 64x32 spatial pix- 
els. We observed in K-band with a spectral resolution of 
FWHM 0.5 nm. The first 12 integrations lasted 5min 
each. During those we noticed NIR activity of Sgr A* 
and we switched to 4 minute exposures. We followed 
Sgr A* for another 32 exposures. In total we interleaved 
nine integrations on a specifically chosen off field (712" 
W, 406" N of Sgr A*). The seeing was wO.5" and the 
optical coherence time ~ 3 ms, some short-time deterio- 
rations excluded. The AO was locked on the closest opti- 
cal guide star (m^ = 14.65, 10.8" E, 18.8" N of Sgr A*), 
yielding a spatial resolution of « 80 mas FWHM, close 
to the diffraction limit of UT4 in K-band (w 60 mas). 

Our detection triggered immediate follow-up obser- 

^This work is based on observations collected at the European 
Southern Observatory, Paranal, Chile. Program ID: 075.B-0547(B) 
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Fig. 1.— The Galactic Center as seen with SINFONI. Chan- 
nel maps from the bright and the dim sample obtained by taking 
the mean in spectral dimension from 2.05 ^m to 2.15 /xm and from 
2.25 /xm to 2.35 fim. All images are scaled the same way and use 
an identical color map. Sgr A* is brighter in the longer wavelength 
maps, indicating that it is is redder than the field stars. Further- 
more, the effect is more pronounced in the dim sample, meaning 
that Sgr A* is redder therein than in the bright sample. 

vations with VISIR, a mid-infrared (MIR) instrument 
mounted at ESO-VLT Melipal (UT3). From 5:25 UT 
onwards VISIR was pointing to the GC. At the position 
of Sgr A* no significant flux was seen. A conservative 
upper limit of 40mJy (not dereddened) at 8.59 /im is re- 
ported (Lagage et al., in prep). 

The reduction of the SINFONI data followed the stan- 
dard steps: From all source data we subtracted the re- 
spective sky frames to correct for instrumental and at- 
mospheric background. We applied flatfielding, bad pixel 
correction, a search for cosmic ray hits, and a correction 
for the optical distortions of SINFONI. We calibrated 
the wavelength dimension with line emission lamps and 
tuned on the atmospheric OH-lines of the raw frames. 
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Fig. 2. — Lightcurve and variation of the spectral power law index 
13 during the flare. Time is counted from 2:40 UT. Top: Flux ratio 
flare/S2. Thin dots are exposures affected by bad seeing. Mid- 
dle/top: /3 using the small apertures background. Middle/bottom: 
13 using the off state subtraction background. Bottom: f3 using the 
constant subtraction background. 

Finally we assembled the data into cubes with a spatial 
grid of 12.5 mas/pix. 

3. ANALYSIS 
3.1. Flux determination 

For all 44 cubes we extracted a collapsed image (me- 
dian in spectral dimension) of a rectangular region 
(0.25" X 0.5") centered on Sgr A* and containing the 
three S-stars S2, S13, and S17. We determined the flux 
of Sgr A* from a fit with five Gaussians to each of these 
images. Four Gaussians with a common width describe 
the four sources. The fifth Gaussian (with a width 3.5 x 
wider, typical for the halo from the imperfect AO cor- 
rection) accounts for the halo of the brightest star S2 
{K w 14 mag). The halos of the weaker sources (all 
K < 15 mag) could be neglected for the flux measure- 
ment. We fixed the positions of all sources (known a- 
priori from a combined cube) and the amplitude ratios 
for the stars. Five parameters were left free: An over- 
all amplitude, the background, the width, the flux ratio 
halo/S2, and the flux ratio Sgr A*/S2. This proce- 
dure disentangles real variability from variations in the 
background, the Strehl ratio, and the seeing. 

As a crosscheck we determined in a second way for 
all images; for both Sgr A* and S2 we measured the flux 
difference between a signal region centered on source and 
a reasonable, symmetric background region. The such 
determined values agreed very well with the fits. For 
further analysis we used the fitted ratios and included 
the difference between the two estimates in the errors. 

3.2. Choice of background 

The value of the spectral power law index f3 crucially 
depends on the background subtraction. Subtracting too 
much light would artificially make the signal look redder 
than it is. Hence, a reasonable choice for the background 
region excludes the nearby sources S2 and S17. Further- 
more, the background flux is varying spatially. Actually 



Sgr A* lies close to a saddle point in the background 
light distribution, caused by S2 and S17 (see Fig. P). 
In the East- West direction the background has a maxi- 
mum close to Sgr A* and in the North-South direction 
a minimum. A proper estimate of the background can 
be achieved in two ways: a) working with small enough, 
symmetric apertures, and b) subtracting from the signal 
an off state obtained at the position of Sgr A* from cubes 
in which no emission is seen. We used the first method 
as well as two variants of the second. 
Small apertures: The local background at a position 
X can be estimated by averaging over a small, symmetric 
region centered on x. Given the background geometry 
we have chosen a ring with inner radius 3 pix and outer 
radius 7 pix. The circular symmetry was only broken 
since we explicitely excluded those pixels with a distance 
to S17 and S2 smaller than 3 and 7 pixels respectively. 
Unfavorable of this method is that the local background 
is only approximated, since a sufficiently large region has 
to be declared as signal region. 

Off state subtraction: The local background can be 
extracted from cubes in which no signal is seen. Since 
the seeing conditions change from cube to cube, one still 
has to correct for the varying amount of stray light in the 
signal region. We estimate this variation by measuring 
the difference spectrum between signal and off cube in a 
stray light region. The latter must not contain any field 
stars and should be as far away from the nearby sources 
as Sgr A*. We used two stray light regions to the left and 
to the right of Sgr A*, between 5 and 10 pixels away. The 
disadvantage of this method is that one needs a suitable 
off state. The latter point is critical for our data. Even 
though Sgr A* has not been detected directly in the flrst 
three cubes, the light at its position appears redder than 
the local background"'^. Assuming that this light is due 
to a very dim, red state of Sgr A*, we would subtract too 
much red light and artificially make the fiare too blue. 
In this sense, this method yields an upper limit for (3. 
Constant subtraction: The off state method can be 
varied to obtain a lower limit for (3. Assuming that the 
true off state spectrum has the color of S2, one can de- 
mand that it is flat after division by S2. In our data, 
the S2 divided off state spectrum is rising towards longer 
wavelengths. Hence, in this third method we estimate 
the background at blue wavelengths and use that con- 
stant as background for all spectral bins. 

3.3. Determination of spectral power law index 

We obtained spectra as the median of all pixels inside 
a disk with radius 3 pix centered on source minus the 
median spectrum of the pixels in the selected off region. 
In order to correct for the interstellar extinction we then 
divided by the S2-spectrum (obtained in the same way 
as the signal). Next the temperature of S2 is corrected 
by multiplying by the value of vSi, of a blackbody with 
T « 25000 K. After binning the data into 60 spectral 
channels it is finally fit with a power law ~ i/^ . With 
this deflnition, red emission has /3 < 0. The error on P 
is obtained as the square sum of the formal fit error and 

^ Inspecting older (non-flare) SINFONI data we found cases sim- 
ilar to the new data and other cases in which the light was identical 
to the local backgro und emission . Thi s is con sistent with the L- 
band observations bv lGhez et"an |20^|20053). 
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Fig. 3. — Correlation plots between the flare flux and the spectral index /3. Points with error bars represent the flare, blue dots are 0si7- 
Thick black dots mark the points with an error A/3 < 1, red triangles mark the data with good seeing ( FWHM < 75 mas) . The line is a 
fit to the thick black dots. Open circles denote the data from lEisenhauer et alj 120051) - near 2 mjy - and lGhez etan 12005bl) - near 7mjy. 
Left: Small apertures method; middle: off state subtraction method; right: constant subtraction method. 

the standard deviation in a sample of 20 estimates for (3 
obtained by varying the on and off region selection. 



4. RESULTS 

We observed a strong (flux density up to 8mJy or 
vL^ K, ergs/s), long (more than 3 hours) flare which 
showed significant brightness variations on timescales as 
short as fO minutes (Fig. El top). While the data is not 
optimal for a periodicity search (poorer sampling than 
our previous imaging data), it is worth noting that the 
highest peak in the periodogram (significance of « 2 cr) 
lies at a period of « fSmin. This is also the timescale 
found by Genzel ct al. (2003), who identify the quasi- 
periodicity with the orbital time close to the last stable 
orbit (LSO) of the MBH. 

We divided the data into three groups: a) the cubes of 
the first peak near t — 50min (" preflare" ) , b) the cubes 
at t > 100 min with T< 0.25 ("dim state"), and c) the 
cubes a.tt> f 00 min with 0.25 ("bright state"). For 
the three sets we created combined cubes in which we 
determined /3, using all three background estimates. In 
all cases we obtained the correct spectral index 2.9 ± 0.5 
for S17 (a star with a spectrum similar to S2 but a flux 
comparable to Sgr A*). For Sgr A* we get: 





preflare 


dim state 


bright state 


off state subtr. 


-1.4 ±0.4 


-0.7 ±0.4 


+0.4 ±0.2 


small apertures 


-1.8 ±0.3 


-1.7±0.4 


-O.f ±0.3 


constant subtr. 


-3.4 ±0.4 


-2.3 ±0.3 


-0.3 ±0.2 



The values obtained from the small apertures lie be- 
tween the values from the other two methods, consistent 
with the idea that they yield upper and lower limits. 
The absolute values vary systematically according to the 
chosen background method. However, independent from 
that, it is clear that the preflare is redder than the dim 
state which in turn is redder than the bright state. An 
obvious question then is whether flux and P are directly 
correlated. 

Hence, we applied the spectral analysis to the individ- 
ual cubes. We kept the data in which Sgr A* is detected, 
the error A/3 < 1.5 and the spectral index for S17 does 
not deviate more than 1.5 from the expected value. The 
resulting spectral indices appear to be correlated with 
the flux (Fig. |3 |3J. The values match the results in 



lEisenhauer eTal] l|2005) and lGhez et all l|2005bj) . Bright 
flares a re indeed bluer than weak flares, as suspected by 
iGhez e t al. ( 2005b) a nd consistent with th e earlier multi- 
band observations of 'Genze l et al.l l)200 3?l. Our key new 
result is that this even holds within a single event. 

For all three background methods it is clear that the 
main event was preceded by a weak, red event. For the 
small apertures and the constant subtraction method in- 
stantaneous flux and spectral index are correlated. For 
the off state subtraction method one could instead group 
the data into a preflare at the beginning and a bluer, 
brighter main event. 

We checked our data for contamination effects. If stray 
light would affect the flare signal, /? should be correlated 
with the seeing (measured by the width from the multiple 
fits). Since we did not find such a correlation, we exclude 
significant contamination. 

5. INTERPRETATION 

Theoretical models predict that the mm-IR emission 
from Sgr A* is synchrotron emission from rclativistic 
electrons cl ose to the LSO (Liu fc Me lia 2001; Quatae^ 
l2003HLiu et al. 2004, 200^. Radiatively inefficient accre- 
tion flow (RIAF) models with a thermal electron popu- 
lation (Te ~ 10 K) produce the observed peak in the 
submm but fail to produce enough flux at 2 fim. The 
NIR emission requires transiently heated or accelerated 
electrons as proposed by Markoff ct al. (2001). 

A conservative interpretation of our data is suggested 
by Fig. 121 middle. There was a weak, red event before 
and possibly independent from the main flare which was 
a much bluer event. Plausibly the preflare i s then due to 
the h igh-energetic tail of the submm peak ijGenzel et alJ 
|2003). The main flare requires nevertheless a population 
of heated electrons. 

A more progressive interpretation follows from the cor- 
relation between flux and /? (Fig. left & right). It 
suggests that the NIR variability is caused by the combi- 
nation of transient heating with subsequent cooling and 
orbital dynamics of relativistic electrons. In the following 
subsections we will exploit this idea. 



5.1. Synchrotron emission 
In the absence of continued energy injection. 



syn- 



chrotron cooling will suppress the high energy part of the 
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electron distribution function. This results in a strong 
cutoff in the NIR spectrum with a cutoff frequency de- 
creasing with time. At a fixed band one expects that 
the light becomes redder as the flux decreases. This can 
cause the correlation between flux and (3. The observed 
flare apparently needs several heating and cooling events. 

The synchrotron cooling time is comparable to the ob- 
served timescale of decaying flanks as in Fig. O (top) 
or IGenz el et al. (2003). In a RIAF model ( Yuan et al.l 
120031) the magnetic field B is related to the accretion 
rate M. For disk models with M w 10"*^ Mg/yr (jAeoll 
I2000t iQuataert fc Gruzinovl 1200(H) one has B « 30 G at 
a radial distance of 3.5 i?s (Schwarzschild radii). The 
cooling time for electrons emitting at A = 2 A2 imtl for 
B = 3OB30G is 

Uooi ~ 8 -830^^^ min . (1) 

This is similar to the orbital timescale, making it difficult 
to disentangle flux variations due to heating and cooling 
from dynamical effects due to the orbital motion. 

5.2. Orbital dynamics 

The timescale of « 20 minutes for the observed vari- 
ations suggests orbital motion close to the LSO as 
one possible cause. Any radiation produced propa- 
gates through strongly curved space-time. Beaming, 
multi ple images, and Doppler shi fts have to be consid- 
ered pollvwood et al.lll997lll999D . Recent progress has 
been made in simulating these effects when observing 
a spatially limited emission region orbiting the MBH 
(jPaumard et a l. 2005: Broderick k L oeb '200 5aiblH) . 

That the emission region is small compared to the ac- 
cretion disk can be deduced from X-ray observations. 



Agol, E. 2000, ApJ, 538, L121 
Bonnet, H. et al. 2004, ESO Messenger, 117, 17 
Broderick, A. & Loeb, A. 2005, MNRAS, 363, 353 
Broderick, A. & Loeb, A., ApJ, 636, L109 

Broderick, A. & Loeb, A., MNRAS, submitted jastro-ph/0509237t 

Eckart, A. et aL 2004, A&A, 427, 1 

Eckart, A. et aL, A&A, in press I astro-ph/0512440"l 

Eisenhauer, F. et aL 2003, Proc. bFlK, 4841, 1548 

Eisenhauer, F. et aL 2005, ApJ, 628, 246 

Genzel, R. et aL 2003, Nature, 425, 934 

Ghez, A. et aL 2003, ApJ, 586, L127 

Ghez, A. et aL 2004, ApJ, 601, L159 

Ghez, A., Salim, S., Hornstein, S., Tanner, A., Lu, J., Morris, M. 

k. Becklin, E. 2005, ApJ, 620, 744 
Ghez, A. et aL 2005, ApJ, 635, in press 
Hollywood, J. & Melia, F. 1997, ApJS, 112, 423 
Hollywood, J., Melia, F., Close, L., McCarthy, D., Dekeyser, T. 

1999, ApJ, 448, 



lEckart et alJ ()2004l |2005() report similar timescales for 
X-ray- and IR-flares. This excludes that the X-ray emis- 
sion is synchrotron light as one would expect cooling 
times < 1 min (eq. ^ . The X-ray emission is natu- 
rally explained by IC scattering of the transiently heated 
electrons that have a 7-factor of 7 « 10"^. The IC 
process scatters the seed photon field up according to 
^ic ~ 7^ fseed- Givcu Vic ~ 10^^ Hz for X-rays the seed 
photon frequency is Vseed ~ 10^^ Hz - the submm regime. 
The X-ray intensity for the synchrotron-self-Compton 
case is given by 

XI IR = Lssc/isync « 10-^^35 W RsV^ B^^^ , (2) 

where ^35 is the seed photon field luminosity in units 
of 10"^^ergs/s and R the radius of the emission region. 
The highest value of L35 compatible with the MIR- and 
submm-data is reached for a flat spectrum (/3 — 0) of 
the heated electrons. Thus, L35 < 1. Via equation (jSJ 
we have R/Rs « O.l^/L^^/^X/IR). From the observed 
ratios X/IR « [O.I - I] we derive R < 0.3 i?s- That 
means, the emission region has to be a small spot. 

A plausible dynamical model is discussed in Paumard 
et al. (in prep). It repr oduces typical lightcurves as ob- 
served bv iGenzel et all (j2003H or in Fig. |21 Due to the 
Doppler effect the observed light corresponds to different 
rest frame frequencies depending on the orbital phase. 
If the source spectrum is curved, ffux and spectral in- 
dex appear correlated. Following this interpretation, the 
emission during the brightest part originates from a rest- 
frequency with larger /? than the dimmer state emitted at 
shorter wavelengths. Such a concavely curved spectrum 
is naturally expected from the synchrotron models. 



Liu, S. & Melia, F. 2001, ApJ, 561, L77 
Liu, S., Petrosian, V. & Melia, F. 2004, ApJ, 611, LlOl 
Liu, S., Melia, F. & Petrosian, V. 2006, ApJ, 636, L798 
Markoff, S., Falcke, H., Yuan, F. & Biermann, P. 2001, A&A, 379, 
L13 

Melia, F., Liu, S. & Coker, R. 2000, ApJ, 545, L117 
Melia, F., Liu, S. & Coker, R. 2001, ApJ, 553, 146 
Narayan, R., Yi, I., & Mahadevan, R. 1995, Nature, 374, 623 
Paumard, T. et al. 2005, AN, 326, 568 
Quataert, E. & Narayan, R. 1999, ApJ, 520, 298 
Quataert, E. & Gruzinov, A. 2000, ApJ, 539, 809 
Quataert, E. 2003, ANS, 324, 435 
Schodel, R. et al. 2002, Nature, 419, 694 
Yuan, F., Quataert, E. & Narayan, R. 2003 ApJ, 598, 301 



